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We develop a sequence of models describing information transmission and decision dynamics for 
a network of individual agents subject to multiple sources of influence. Our general framework is 
set in the context of an impending natural disaster, where individuals, represented by nodes on 
the network, must decide whether or not to evacuate. Sources of influence include a one-to-many 
externally driven global broadcast as well as pairwise interactions, across links in the network, in 
which agents transmit either continuous opinions or binary actions. We consider both uniform and 
variable threshold rules on the individual opinion as baseline models for decision-making. Our results 
indicate that 1) social networks lead to clustering and cohesive action among individuals, 2) binary 
information introduces high temporal variability and stagnation, and 3) information transmission 
over the network can either facilitate or hinder action adoption, depending on the influence of the 
global broadcast relative to the social network. Our framework highlights the essential role of local 
interactions between agents in predicting collective behavior of the population as a whole. 



INTRODUCTION 

The influences that humans have on one another's 
opinions and subsequent actions can shape large-scale 
movements [US], and are often facilitated by the shar- 
ing of information. Advances in information technologies 
are rapidly changing the way that humans exchange and 
share information. The widespread adoption of radio in 
the 1920s and television in the 1950s ushered an era of 
live broadcast media, greatly increasing the speed and 
scope of mass communication from that of newspapers, 
which had been the dominant form of "one-to-many" 
information dissemination since the development of the 
Gutenberg printing press in the 15th Century. The com- 
mercialization of the Internet in the late 1990s catalyzed 
the development of new digital broadcast services, con- 
sumed by a diversity of computing machines that now 
include laptop computers, mobile "smart" phones, and 
other handheld devices. These devices are now ubiqui- 
tous; by mid 2010 there were more than 5 billion mobile 
phone connections worldwide [5], with some regions ex- 
periencing more than 100% penetration (meaning that 
there is more than one mobile device per person). 

A distinguishing feature of the Internet and these mod- 
ern digital devices is that they enable rapid "one-to- 
many" and "many-to-many" communication, using ser- 
vices such as Facebook and Twitter, a phenomenon that 
has come to be known as social media. By the end of 
2011 more than 300 million users were accessing Face- 
book using mobile devices [6j [7]. In this modern era, 
anyone with a digital device and an Internet connection 
can publish information, and this is dramatically chang- 
ing the roles that individuals and corporations play in 
traditional media industries such as books, music, film, 
and news journalism [8]. 

There is a general recognition that these technologies 



allow information to spread faster and perhaps more ef- 
fectively. Such social epidemics, like biological epidemics 
[9j [10] , can be thought of as the result of single- [TT] or 
multi-stage [12] complex contagion processes [T3HT5] that 
can propagate through Facebook [16] , news websites [T7] 
and other social media [18] . blogs [19] and Twitter [20] . 

Recent events on the East Coast of the United States 
have shown that social media can be helpful during ex- 
treme weather events and natural disasters [21]. In the 
case of the 5.8 magnitude earthquake that occurred in 
Virginia on 23 August 2011, the social network service 
Twitter proved to be more reliable than cellular phones, 
which became overloaded by increased call volume imme- 
diately after the event [22], and news of the event prop- 
agated to nearby New York on Twitter faster than the 
seismic waves themselves [23 . Social media also allows 
for rapid organization of disaster response information. 
For example, during the 7.1 magnitude earthquake in the 
Canterbury region in the South Island of New Zealand 
on 4 September 2010, the University of Canterbury used 
Facebook and other web-based technologies as a promi- 
nent source of up to date information and support for 
many months [24]. It is less clear how this faster infor- 
mation propagation affects the collective behavior of the 
individuals who consume it. Social networking technolo- 
gies have been credited with facilitating the 2009 revo- 
lution in Moldova [25] and the "Arab Spring" uprisings 
in 2011 [26, 27 . More prominently, social media has be- 
come an important component of corporate marketing 
and advertising, with considerable effort now directed at 
determining how to optimize the use of this new media 
[28H30]. 

In this article, we consider the following question: Do 
social networking technologies like Facebook and Twitter 
"help" to bring a group of individuals to action? We con- 
sider the specific case of a population of individuals, each 
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of whom must decide if and when to commit to a binary 
decision, and we assume that these individuals are ex- 
posed to information both from broadcast sources and 
over social media. The study of information diffusion 
on social networks has a lengthy history of illuminating 
the large-scale spread of rumors, social norms, opinions, 
fads, and beliefs [9j [13j I3TH35] . The particular set of 
problems where individual decision-making occurs in the 
presence of external information is known as "binary de- 
cisions with externalities" [3TJ |36j [37] . Here the individ- 
ual's decision is often modeled as a threshold on their un- 
derlying opinion [3TJ [38ti4Tj , which is modulated by the 
opinions of other individuals. We pose a minimal model 
of individual decision-making and consider the effect of 
different forms of information exchange on the behavior 
of the group as a whole. 

As a concrete example, we focus on a situation where a 
population is faced with a pending natural disaster (e.g., 
hurricane or wildfire). In such circumstances, it is know 
that each individual culls information from a wide vari- 
ety of digital, sensorimotor, and social sources [42] and 
must decide if and when to evacuate. We assume that 
each individual has a simple decision rule: if the individ- 
ual believes that the disaster is sufficiently likely, then 
he or she will evacuate. Individuals receive information 
about the disaster from a "global source" that broad- 
casts updates to the population as a whole (see Fig.[l]A.), 
and these individuals also exchange information over a 
social network that allows them to share opinions and 
observe the binary decisions of others (i.e., who has evac- 
uated; see Fig[l|3). By numerically exercising a series of 
increasingly complex models, we illustrate the tensions 
and trade-offs inherent in social decision dynamics. 

Our results suggest that information transmission over 
the social network can either facilitate or hinder the ac- 
tion adoption depending on the influence of the global 
source relative to the social network. Further, we find 
that the sharing of binary information results in high 
variability, cascade-like dynamics in which the time of 
collective action is difficult to predict. 



MODEL CONSTRUCT 

We build a model of decision-making and social net- 
work interaction in discrete time. The social networks we 
study can be thought of as Facebook- or Twitter-like, in 
the sense that an agent posts updates or tweets to these 
media and subsequently other agents can check the post- 
ing at some rate. In our model framework, each agent 
receives a directed update from the other agents one at 
a time at prescribed rates. 

Let t = 0, 1, 2, . . . index the discrete time increments. 
The social network consists of n individuals, in which 
each agent j = 1, 2, . . . , n has two state variables: 



• Sj(t) = the internal state at time £, where Sj(t) G 

[o,i]; 

• Xj(t) = the externally observable state at time t, 
where Xj(t) G {0,1}. 

The internal state assumes continuous values, but the ex- 
ternally observable state Xj(t) reflects a binary decision 
on the part of agent j, derived from the decision rule 
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if Sj(t) > tj 

HS j (t)<T j 



(1) 



and where Tj G [0, 1] represents a threshold value for 
agent j. The use of a threshold for the binary decision 
Xj(t) is consistent with long-standing modeling efforts 
of collective behavior beginning with the work of Gra- 
novetter in the 1970s [38] and is particularly relevant for 
decisions which are inherently costly [31] (such as evacu- 
ations). 

We think of the internal state information as "private" 
in the sense that, for example, it reflects an underlying 
belief on the part of the individual. While it is possible 
that an individual could share this private information 
with another, this exchange requires the individual to 
volunteer it. In contrast, the externally observable state 
is "public" in the sense that anyone who sees that indi- 
vidual would observe their binary decision. It has been 
argued that this separation between internal opinion and 
external decision is critical for an understanding of how 
convictions are coupled to actions [43[ [44] . 

We represent the internal state of all agents at time t 
using S(t) = {Si(t), 52 (t), . . . , 5 n (t)}, and we represent 
the externally observable state of all agents at time t 
using X(t) = {X 1 (t),X 2 (t),...,X n (t)}. 

The internal state of an agent changes over time as 
a result of information from three different sources: 1) 
a global source via a "global broadcast", 2) the inter- 
nal state information of friends via "social sharing" , and 
3) the decision state of neighbors via "neighbor observa- 
tion". Each of these sources is described in more detail 
below. 

Global Broadcast. We introduce a special external 
agent called the "global source" of information that in- 
fluences, but is not influenced by, the other agents. Let 
G(t) G [0, 1] represent the value that is broadcast by the 
global source in time period t. We assume that receipt 
of a broadcast message by agent j is binary (i.e., either 
it happens or it does not). Let Uj(t) = 1 represent the 
successful broadcast from the global source to agent j 
at time t, and Uj(t) = represent no broadcast. Thus 
the vector U(t) = {ui(t) , v,2(t) , ...,u n (t)} represents the 
overall broadcast from the global source at time t. This 
global source is the primary "external driver" of dynam- 
ics in this system. 

Social Sharing. We assume that the action of social 
sharing between agents is binary, either occurring or not 
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occurring at any point in time. That is, we let a^-(t) = 1 
if agent i shares its internal state information with agent 
j in time period t, and = if not. Thus, the matrix 

A(t) = {a>ij(t)} represents the adjacencies for exchange 
of internal state information among the agents in time 
period t. By construction, CLjj(t) = for all agents j and 
all time periods t. In general, we use symmetric shar- 
ing (i.e., dij(t) = ciji(t)) to indicate that any agent may 
address any other agent, however the mathematics here 
do not require it. In our model, however, these links are 
directed in the sense that for a given sharing event, an 
agent asks for another agent's opinion or state, but does 
not share its own state. Whether or not there is sharing 
symmetry or directionality depends on the specific ap- 
plication we hope to model. Communication across cell 
phones tends to be symmetric and undirected while that 
across Twitter and Facebook are not necessarily so. 

Neighbor Observation. Similarly, we assume that the 
observation of another agent's external state either hap- 
pens or does not in each time period t. Let B(t) = 
{bij(t)} represent an adjacency matrix for exchange of 
externally observable binary state information among the 
agents, where bij G {0,1}. By construction, bjj(t) = 
for all agents j and all time periods t. The matrix B 
might represent a network of interactions based on phys- 
ical location, where nearby individuals literally see one 
another even if they do not communicate. For purposes 
of our discussion, we refer to agents who are adjacent 
in the matrix B as neighbors. That is, we narrowly de- 
fine two agents as "friends" if they share internal state 
information, and we define two agents as "neighbors" if 
they can observe the external state of each other. Again 
the mathematical formulation here can be used to model 
both symmetric and non-symmetric sharing. 

General Update Rule. In the presence of all three 
sources of information, the internal state of agent j 
evolves according to the general update rule 

Sj(t+1) = 

$>ii(*)$(*) +Y,b ij (t)X i {t) + Uj(t)G(t) 
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FIG. 1. Model Construct Here we illustrate our multi- 
layered broadcast-social network system composed of a graph 
allowing information diffusion between a global source and a 
population of 12 individuals (broadcast network) and a sec- 
ond graph allowing information diffusion between individuals 
(social network). (A) The broadcast network structure is a 
directed bipartite graph where the link weights are equal to 
the probabilities that the global source transmits information 
to the individual (probabilities given by the vector p, which 
varies over multiple numerical simulations M). (B) The so- 
cial network structure is an all-to-all undirected graph where 
the link weights are equal to the probabilities that an individ- 
ual transmits information to another individual (probabilities 
given by the matrix P, which varies over multiple numeri- 
cal simulations). In concrete terms, Pij can be thought of as 
the combined rate that agent i posts and agent j reads the 
posting on Facebook. 



While opinion dynamics have been studied using a va- 
riety of models [45], the averaging rule given in Eq. (§ 
is most consistent with the Hegselmann-Krause model of 
opinion dynamics [46 , models of "continuous opinion dy- 
namics" [47], and the many models of coordination and 
consensus of autonomous agents [37 , 48 - 51. The discrete 
nature of this update rule is consistent with the fact that 
information is often issued at some frequency or can be 
obtained in discrete units from governmental, social, or 
technical sources |42l [52] . 
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This update rule is a deterministic averaging of the cur- 
rent internal state of agent j, the internal state of agent 
j's friends, the external state of agent j's neighbors, and 
any global broadcast information. However, we assume 
that the coefficients in A, £?, and U are stochastic in 
time; that is, in any given time period t the information 
to agent j might or might not be received. This stochas- 
ticity is consistent with the fact that communication over 
these types of networks is likely to be sparse - due to geo- 
graphic and energetic constraints - and dynamic - due to 
agent movement and constraints on communication [37] . 



COLLECTIVE DECISION DYNAMICS 

In what follows, we are interested in the real-time ex- 
change and mixing of information among the individuals, 
and how it affects the decision-making of the collective 
group. 

We measure the behavior of the system as a whole in 
several ways. Recall that the social network consists of n 
individuals, in which each agent j = 1, 2, . . . , n has two 
state variables: Sj (t) and Xj (t). We refer to the sequence 
S(i),t > which indexes over the j agents as the infor- 
mation state trajectory for the system, and the sequence 
X(t),t > similarly as the action adoption trajectory. 
Let m = 1, 2, . . . , M index the numerical trials associated 
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with a particular experiment, and let £™(£) represent the 
value of Sj (t) during the m th trial. We compute the aver- 
age information state of the population at time t during 
experiment m as (*S™(£)) n = ^ Ej=i SJ 1 ^)] accordingly, 
the sequence (S™(t)} n ,t > is the average information 
state trajectory for the system during experiment m. We 
compute the average information state of individual j at 
time t as {Sf(t)) m = ^Eli^W; this is the aver- 
age information state of any individual in the population. 
Accordingly, the sequence (S™(£)) m ,£ > is the ensem- 
ble information state trajectory of a single individual over 
the ensemble of numerical trials. Finally, we can estimate 
the average ensemble information state trajectory as 



given by 



m ( n \ 
m=l \ j=l J 

j — l \ rn—\ / 



(3) 



It is important to note that the variance of Sj(t) over 
individuals could be very different from the variance ex- 
pected over numerical trials; for example, one distribu- 
tion might be normal and the other heavy-tailed. The 
separation of the individual and ensemble averages allows 
us to separately probe these distinct sources of variability 
in the system. 

We introduce the term N(t) = Y^j=\Xj(t), which 
counts the number of individuals whose internal state 
has exceeded their decision threshold and who therefore 
have taken action. With this definition, we observe that 
(Xj(t)} = N(t)/n measures the fraction of individuals 
who have chosen the binary action. 

In many cases, we are also interested in the amount 
of time it takes for the population to collectively adopt 
new information or take action. Define Hs(ot) = 
min{t | (Sj(t)) > a} to be the first hitting time for the 
population to reach some average information level a £ 
[0, 1]. Similarly, define Hx(ct) = min{£ | (Xj(t)) > a} to 
be the first hitting time for the population to reach some 
average adoption level a £ [0,1]. For example, Hx(0.5) 
is the amount of time for half the population to adopt 
the binary action. 

Conceptually, we can think of two basic types of in- 
formation diffusion; see Fig. [I] The first is the (one-to- 
many) broadcast of external information by the global 
source to the agents; we can think of this interaction 
in terms of a broadcast network. The second (many-to- 
many) diffusion occurs on two different social networks: 
one network of friends and another network of neighbors. 
We study the behavior of the system in the following spe- 
cific cases. 

Case 1: Information broadcast from the global source 
only. In this case, the only information exchange is the 
broadcast from the global source, and the update rule is 



Sj(t- 



1) 



S^ + u^Gjt) 

1 + Uj(t) 



(4) 



In the simple case where Sj (0) = and where the global 
source broadcasts G(t) = 1 for all time periods t, then 
the progression of internal state for agent j goes as 1 — 
l/(fc + l) after the k th update. That is, with each update 
from the global source, agent j moves "half way" to the 
broadcast value of 1, reaching it only in the limit as t — > 
oo. 

As noted above, we assume that broadcast messages 
from the global source are received stochastically. This 
could occur in practice because agent j only "tunes in" to 
the global source sometimes or because the global source 
has limited success in its ability to reach agent j. For 
each time period t we generate a vector U (t) as an inde- 
pendent and identically distributed random vector, where 
Prob{uj = 1} = Pj, and this probability is independent 
for each agent j. Thus, E[iZj(t)] = Pj captures the ex- 
pected broadcast "rate" [52] to agent j (assumed to be 
stationary for now), and E[?7 (£)] represents the overall ex- 
pected broadcast from the global source in time period 
t. 

In this simple case, we can derive expected values for 
Sj(t) and Xj(t) analytically. Specifically, after t discrete 
time units, the probability of agent j having received 
k < t updates is given by the binomial distribution, 
Binom(fc; The expected value of Sj(t) is therefore 



k=0 v 7 \ 1 / 



(5) 



In addition, E[Xj(t)] is simply the probability that agent 
j has taken action (i.e., has adopted state 1) and can be 
defined as 



E[X j (t)]=Prob{S j (t)>T j } 



(6) 



where I| a } is the indicator function, namely I| a | = 1 
when condition a is true, and I{ a j = when condition a 
is false. 

The values E[Sj(t)] and E[X,-(t)] in the simple case 
of global broadcast serve as a baseline against which we 
can evaluate the impact of various types of social network 
exchange. 

Case 2: Global broadcast with social network exchange. 

We now add information exchange among friends 
alongside the global broadcast. This gives us the fol- 
lowing update rule: 



Sj(t+1) 



^2a zj m(t)^ Uj (t)G(t) 

i 



(7) 
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Again, we assume that exchange of information be- 
tween friends is stochastic in time. Thus, we generate 
each matrix A(t) as a weighted matrix, where Prob{aij = 
1} = Pij and this probability is independent for each 
pair. We let this probability be stationary in time, 
based on the common assumption that the information 
transmission occurs on a faster time scale than network 
changes [53]. Thus, E[a^-(t)] = P^ represents the ex- 
pected rate at which agent i influences agent j, and 
E[A(t)] is the expected information exchange within the 
social network of individuals, as might happen using tech- 
nologies such as mobile phones, Facebook, or Twitter. 

This case corresponds to a particular type of consensus 
problem for which there are analytic results that describe 
the convergence of E[Sj(t)], specifically the conditions 
under which it is guaranteed and how long it will take. 
Jadbabaie et al. [49] consider the case of "leader follow- 
ing" in consensus problems in which one of the agents 
never updates its own variable, but indirectly influences 
all of the other agents. The role of the "leader" is equiv- 
alent to that of our global source. (Bertsekas and Tsit- 
siklis [54 argue that this result is essentially a special 
case of the more general result in [55 .) Similarly, Jad- 
babaie [56] discusses routing in networks in which nodes 
iteratively update their coordinate information, but cer- 
tain "boundary" nodes retain fixed locations (again, like 
our global source). A third related example is the model 
of Khan et al. [57 , which describes consensus on random 
networks in which there are "anchor nodes" (again, like 
our global source in that their state does not change) and 
"sensor nodes" (like our agents who update through ran- 
dom mixing). Finally, Galam et al. [58, 59 use the term 
"inflexible agents" to indicate nodes whose opinions do 
not change and find in some cases that these agents drive 
opinion dynamics. 

Case 3: Global broadcast with binary information only 
among neighbors. In this final case, we consider the im- 
pact of binary exchange among neighbors alongside the 
global broadcast. This gives us the following update rule: 

+ + «*(*) 

Sj(t+1) = ^ , (8) 

1 + 6«(t) + «,•(*) 

i 

Again, we assume that this type of exchange is stochas- 
tic, and we generate each matrix B(t) as a weighted ma- 
trix, where Prob{bij = 1} = Q^-, and this probability 
is stationary and independent for each pair. Thus, 
Effect)] = Qij, and E is the expected information 
exchange within local neighborhoods. In this case, we an- 
alyze the network of friends separately from the network 
of neighbors. In future simulations, it may be interesting 
to examine the special case in which the two networks 
are either identical or change in a dependent manner. 

The behavior of this system is sufficiently complicated 
that analytic results do not, to our knowledge, exist for 



this case. We therefore turn to numerical simulations in 
order to analyze and compare this case to the simpler 
ones described above. 



NUMERICAL EXPERIMENTS 

We conduct several numerical experiments in which 
we simulate the behavior of n = 100 agents over a max- 
imum of T = 1000 time units. Each agent j starts with 
Sj(0) = 0, and the global source broadcasts G(t) = 1 for 
all time periods t. For each of our three cases, we simu- 
late a total of M = 100 trials, and for each trial we select 
fixed broadcast rates for each agent that are uniformly 
distributed on the interval [0,1] (that is, pj ~ U[0, 1]). 

Similarly, for each trial we select social sharing rates 
Pij or neighbor observation rates Qij, each of which re- 
main fixed within an individual trial. We define a sym- 
metric matrix R whose diagonals decrease in value: 

Ri,i+\k\ =(n-k + l)/n V fc e {1, 2, ...(n - 1)}. (9) 

The matrix R can be thought of as a probabilistic form 
of a regular lattice network. We choose P and Q as 
randomly scrambled versions of R, where the scrambling 
maintains the symmetry of the matrix. In general, this 
approach provides flexibility for the examination of topo- 
logical structures between random and regular graphs. In 
this work, we focus on random topologies to understand 
benchmark behavior. The distribution of probabilities 
across agents provides a weighted counterpart to degree 
heterogeneities examined in other studies of opinion dy- 
namics [60j [61] , which have important consequences for 
collective action. 

For each of our three cases, we investigate the evolution 
of information states and action adoption for individual 
agents and the collective group as a whole. 

Average Information State 

Figure [2] compares the evolution of information states 
across the three case studies. Figure [2]A. shows the step- 
wise trajectory Sj(t) for a single agent, and Figure [2^3 
shows the overall population of agents who independently 
update their state information in the presence of the 
global broadcast only. Figure [2p presents the average 
information state trajectory for the collective group cal- 
culated both through simulation (black) and through the 
analytic solution given in Eq. ([5|, which agree well. 

Figures [2p-F present the equivalent figures for the case 
of global broadcast and social network sharing. We ob- 
serve that in the presence of social networking, individual 
agents (Figure [2jZ>) as well as the entire group of agents 
evolve their state information in a collectively smooth 
manner (Figure |2jil). However, the overall rate of increase 
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Case Study 1 : Global Broadcast Only 
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Case Study 3: Global Broadcast with Neighbor Observation 
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FIG. 2. Information Diffusion Under Three Types of Communication. Panels A, D, G: Information state trajectory 
for a single randomly selected agent j as a function of time. Panels B, E, H: Information state trajectories for a population 
of n — 100 agents. Information states (color) for each individual (row) are shown as a function of time (x-axis). Panels C, F, 
I: Average information state variable, (Sj), averaged over M=100 numerical simulations as a function of time. Note that the 
simulated curve (black) in Panel C is accompanied by a curve that was analytically calculated from Eq[5](red). 



for the average information state in the case of social net- 
work sharing (Figure [2^) is slower than in the case of 
global broadcast only (Figure |2j3). 

Figures [2p-I present the equivalent results in the case 
of neighbor observation in which the only agent-to- agent 
information exchange comes from observing the external 
state of other agents. For the given parameters (n=100, 
d=l, M = 100; for a description of d see later sections), 
we observe that none of the individuals ever raises their 
information state significantly from the starting value 
of zero. This case represents an extreme of the previ- 
ous one in which social network sharing through neigh- 
bor observation slows any rising of the information state. 
(Note: See later sections for other areas of the state space 
in which potentially more interesting behaviors like cas- 



cades and stagnation occur.) 



Average Adoption State 

Figure [3] displays the adoption trajectories for the three 
cases, when each individual agent j has a decision thresh- 
old Tj = 0.5. In Figure |3j^, we observe that under 
only global broadcast the information state of individual 
agents evolves in a manner consistent with Figure [2^3, ex- 
cept that we terminate each trajectory when the individ- 
ual crosses the decision threshold. Figure [3^3 shows the 
average adoption state trajectory for a single trial calcu- 
lated both through simulation and through the analytical 
solution given in Eq. ([6]), which agree well. Importantly, 
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to remain comparable to our simulations in which p is 
chosen ~ ?7[0, 1], Eq. ([6| was calculated for a range of 
< pj < 1, and the average is plotted in Figure [3^3. 
We note that the action adoption curve for a population 
with fixed pj is not equivalent to that for a population 
with distributed pj . Finally, Figure |3p shows the effect 
of different threshold values r on the average adoption 
state. 

Figures [2J3-F present the equivalent results for the case 
of global broadcast and social network sharing. As be- 
fore, we observe that in the presence of social networking, 
the entire group of agents evolves their state information 
in a collectively smooth manner (Figure J2p). Since this 
groups shares a common decision threshold, we observe 
that action adoption now occurs abruptly (Figure [2^1) 
and that the onset of this abrupt change depends pri- 
marily on the threshold level r (Figure |2^). 

In the case of neighbor observation and global broad- 
cast, Figures [3]3-I again show that none of the individuals 
ever raises their information state significantly from the 
starting value of zero. These results are consistent with 
those shown in Figure [2p-I and represent an extreme of 
social inertia keeping the information state at very low 
values. (Again, see later sections for other areas of the 
state space in which potentially more interesting behav- 
iors like cascades and stagnation occur.) 

Heterogeneous Decision Thresholds 

Up to this point, we have assumed that individual 
agents share a common decision threshold Tj. In prac- 
tice, this is unlikely to be the case and recent theoretic 
and experimental work suggests individual variability in 
decision thresholds may significantly alter the dynamics 
of the population [62j [63] . We now consider the implica- 
tions of a heterogeneous population of agents where the 
threshold for agent j is uniformly distributed tj ~ J7[0, 1]. 
Figure [4] displays the impact of variation in underlying 
Tj on the mean adoption state of the collective group as 
a function of time when averaged over multiple trials, 
given by (Eq. [3|, as well as the corresponding variance 
of the adoption state of the collective group as a function 
of time. 

We note that in the case of the global broadcast only, 
the distribution of thresholds, whether homogeneous or 
heterogeneous, has little effect on the mean or variance 
of the adoption state trajectories (Figure |I]A-B). 

However, in the case of global broadcast and social 
network sharing, having a homogeneous distribution of 
thresholds over the population leads to a 'tipping point', 
where at one time point no one has taken action while 
a few time steps later, everyone has taken action (Fig- 
ure |4p, black line). For heterogeneous thresholds, this 
drastic adoption is no longer evident, and instead the 
transition from zero adoptions to all adoptions is smooth 



and gradual (Figure [4p, red line). Consistent with these 
results, we find that the variance in the adoption states 
across individuals is large precisely at the tipping point 
for homogeneous thresholds (Figure [4)3, black line) and 
is small for all times when heterogeneous thresholds are 
used (Figure |4jZ), red line), suggestive of the formation of 
a 'collective'. 

In the case of neighbor observation and global broad- 
cast, Figure [4^1 shows that for this set of parameters 
fa = 0.5 OR Tj ~ U[0, 1]), when the threshold is 
distributed homogeneously the population is unable to 
take action (black line) but when it is distributed het- 
erogeneously, the population can take action (red line). 
The fact that the population can not take action when 
Tj = 0.5 is true for all simulations, leading to a zero- 
variance over simulations (Figure [4^, black line). In con- 
trast, the action adoption trajectories for the population 
vary over simulations when Tj ~ U[0, 1], and they do so 
particularly for time points close to the point of maxi- 
mum adoption (Figure [4^, red line). 

Stagnation and Cascades 

The adoption state trajectories reported in the previ- 
ous section were averaged over M numerical simulations. 
However, in order to fully understand the case of neigh- 
bor observation and global broadcast, it is important to 
examine individual simulation trials in addition to their 
average. 

In Figure [5^., we show information state trajectories 
for all individuals in a population in four different sim- 
ulated trials. These four examples highlight the differ- 
ent possible group behaviors including stagnation (Fig- 
ure J5|A1), partial cascades occurring either early (Fig- 
ure |5pV2) or late (Figure [5]A.3) , and full cascades (Figure 
|5]A.4). This high variability is an important factor that 
sets Case 3 (neighbor observation and global broadcast) 
apart from the other two cases. Cascade behavior has 
previously been demonstrated in wide variety of topo- 
logical structures and update rules [31] [4TJ [64-66 and 
understanding when and how they occur if of particular 
import for predicting large-scale social movements. 

We can study this behavior more systematically by ex- 
amining the first hitting time for the population to reach 
an average adoption level of 0.5 (e.g., Hx(0.5)). We de- 
fine a 'converging simulation' as one for which Hx(0.5) 
is identified within 1000 time steps and a nonconverging 
simulation as one for which it is not. In Figure [5j3, we 
show the complementary cumulative distribution func- 
tions (CCDFs) of converging simulations for this model 
(red line). We note that the majority of converging sim- 
ulations reach the target average adoption level in rela- 
tively short times. However, the distribution is heavy- 
tailed such that much longer times are consistent with 
the statistics. In fact, out of a total of 5000 simulations, 
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FIG. 3. Action Adoption Under Three Types of Communication. Panels A, D, G: Adoption state trajectories for 
a population of n — 100 agents. Adoption states (color) for each individual (row) are shown as a function of time (x-axis). 
Note that the color white indicates that the individual has taken action. Panels B, E, H: Average adoption state variable, 
(Xj), averaged over M=100 numerical simulations as a function of time. Panels C, F, I: Dependence of average adoption state 
variable, (Xj), on the threshold level r, which we have varied in this figure from 0.1 to 0.9. Note that the simulated curve 
(black) in Panel B is accompanied by a curve that was analytically calculated from Eq. (|6l (red). 



only 2127 solutions converged in 1000 time steps indi- 
cating that stagnating periods can extend beyond the 
studied temporal window (0 < t < 1000). In the next 
section, we will further explore this behavior as a func- 
tion of the amount of exchange between social and global 
information sources. 



In this section, we introduce an influence parameter 
< d < 1, which serves to tune the network effects in 
Case 2 (social sharing) and Case 3 (neighbor observa- 
tion), respectively, as follows: 



Impact of exchange rates 



The behavior of the system under all three cases high- 
lights the role of different types of information exchange. 
We observe that our stylized form of social network ex- 
change tends to move the group as a whole. 



Sj (t)+d^2 a ij W s i ( f ) + U 3 CO 

Sj(t + M) = ^ (10) 

l + d2^aij(t) + Uj(t) 
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FIG. 4. Action Adoption Under Heterogeneous 
Thresholds. Thresholds are either chosen to remain con- 
stant over all individuals ('homogeneous'; black lines) or to 
vary over individuals ('heterogeneous'; red lines). Panels A, 
C, E: Average adoption state variable, (Xj), averaged over 
M=100 numerical simulations as a function of time. Panels 
B, D, F: The variance (over M — 100 numerical simulations) 
of the average adoption state variable, (Xj), as a function of 
time. Note that the curves shown in panel (D) are on such 
different scales that we have plotted them in separate subplots 
to enhance visualization. 



and 



Sj(t+l,d) = 



l+d^bijW+ujit) 



. (11) 
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FIG. 5. Variability in Action Adoption Behavior (A) 

Example plot of the information state (color) trajectory of 
agents in the population (y-axis) as a function of time (x- 
axis) showing (1) stagnation, (2) partial cascade and then 
stagnation, (3) late partial cascade, and (4) full cascade. (B) 
The complementary cumulative distribution function (CCDF) 
of the first hitting time for the population to reach the average 
adoption level of 0.5 (e.g., if* (0.5)) for a set of M = 5000 
simulations as a function of the influence parameter d. 



network has full influence, but when d = these cases 
each reduce to global broadcast only where the social net- 
work has no influence. This is illustrated in Figure [6]A. 
and F, where we see that for low values of <i, the informa- 
tion state trajectories are similar to those found in Case 1 
(global broadcast only; compare to Figure [3^). The con- 
cept of network influence is akin to the concept of infor- 
mation credibility, which varies over modern information 
transmission technologies [67, 68J: a network with high 
credibility could be modeled as one with greater influence 
and visa versa. 



When d = 1 these cases remain unchanged and the social 



As d increases, we find that the variance in informa- 
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FIG. 6. Effect of the Influence Parameter in Case 2 (Global Broadcast with Social Network Sharing; A-E) and Case 
3 (Global Broadcast with Neighbor Observation; F-J). (A,F) Information states for a population of n=100 individuals as a 
function of time for four values of the influence parameter: box outlines correspond to d- values given in the legends of panels 
(B) and (G). As the influence parameter increases, the variance of information states across individuals decreases (B, G), 
indicating the formation of a collective, and the mean of information states across individuals decreases (C, H), indicating that 
it takes more time for individuals to reach higher values of Sj. Importantly, the four lines shown in (B-C,G-H) are average 
(over M = 100 simulations) trajectories drawn from a continuous measure of influence, d. In (D-E, I-J) we therefore show the 
full average contour plots for 7400 numerical simulations. Four colored bars on the left of each plot indicate the d- values for 
which the four lines in (B-C,G-H) are taken. In these contour plots, color indicates the variance (top) or mean (bottom) of the 
information state trajectory as a function of time (x-axis) for influence parameters in the range 0.0001 < d < 1 (y-axis). 



tion state trajectories across agents decreases (Figures 
[6^3 and G), which is consistent with the formation of 
a single 'collective' state as the social network becomes 
stronger. In addition, as d increases we also find that 
it takes longer for agents to reach any given state value, 
indicating that the social network is maintaining some 
inertia and holding agents closer and closer to their orig- 



inal states (Figures [6p and H). These general behaviors 
are also evident in a more continuous phase diagram of 
the state space of the system (Figures [6p-E and I-J). 
We note that the rightward tail of nontrivial behavior in 
Figure |6)[ in comparison to Figure [6p, and the extension 
of low (Sj) to longer times in Figure [6]j in comparison 
to Figure [6^1 are a result of the stagnation followed by 
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cascading behavior present for the case of neighbor ob- 
servation. These results indicate a wider range of active 
change or variability in model behavior relative to the 
case of social network sharing. 

Importantly, the high variability in population behav- 
ior in Case 3 in the previous section is present over a 
wide range of influence parameter values (Figure [5^3). 
Here we study the first hitting time for the population 
to reach an average adoption level of 0.5 (e.g., i^x(0.5), 
when 50 out of a possible 100 agents have adopted). Ac- 
cumulating data over 5000 numerical simulations, we find 
that it displays heavy tailed behavior: in the majority of 
simulations, it takes a short time for those 50 agents to 
take action, but in a few simulations, it takes a very long 
time. Furthermore, this heavy-tailed behavior is sensitive 
to the influence parameter d and therefore the strength 
of the social network. For low values of more numer- 
ical simulations have shorter first hitting times which is 
consistent with what we would expect for situations that 
approximate Case 1 (global broadcast only). This result 
is consistent with the fact that the social network retards 
information state progress. 

Does the social network help or hinder? 

An important tradeoff is evident in Figure |6p and [6)H[ 
for low values of d. When d = 0.0001, the average in- 
formation state trajectory curve rises the fastest initially 
but then slows for later values of t. In fact, for a slightly 
larger value of d = 0.001, the average information state 
trajectory curve lags the d = 0.0001 initially and then 
surpasses it later. This result suggests that the social 
network - if weakly present - can help the entire popu- 
lation take action sooner. However, when the social net- 
work becomes stronger (e.g., d = 0.01 and d = 1), it acts 
like a cage, keeping information states of all individuals 
from rising swiftly. 

We investigate this behavior more systematically in 
Figure [7| where we observe three distinct behaviors: 1) 
facilitation of information propagation for small rates of 
social exchange (e.g., small values of the influence param- 
eter d] Figure [7^, bottom), 2) hinderance of information 
propagation for large rates of social exchange (e.g., large 
values of <i; Figure [7^, top), and 3) a combination of first 
hinderance and then facilitation for intermediate rates of 
social exchange (e.g., intermediate values of d\ Figure[7^, 
middle). Transitions between these three distinct behav- 
iors are smooth, as demonstrated in Figure [7^. 

These results demonstrate that the social network can 
both hinder and facilitate information state changes and 
by extension action adoption. When the influence of the 
social network is large, the inertia of popular opinion 
dampens the effect of the global source attempting to 
inject new information into the system. On the other 
hand, when the influence of the social network is small, 



the added mixing facilitates the dissemination of the new 
information provided by the global source. 



DISCUSSION 

Our long-term objective is to develop a framework that 
enables predictive modeling of collective decision dynam- 
ics in situations involving multiple sources of information. 
Modern communication technologies and social network- 
ing applications provide fast, global means of informa- 
tion dissemination. The need to determine the impacts 
of these technologies on individual decisions and by ex- 
tension collective action make it essential to understand 
the interplay between multiple information sources. 

This paper lays the foundation of such a framework, 
by systematically exploring a sequence of models that 
aim to capture tradeoffs and tensions that arise when a 
global broadcast source competes with information trans- 
mission between individual agents. Despite our necessar- 
ily simplified scenario, we find that information trans- 
mission over the network can either facilitate or hinder 
action adoption, depending on the relative influence of 
the global and social information sources. In most situa- 
tions the social network acts overall as a damping force, 
homogenizing opinion states and delaying action adop- 
tion in the population. 



Insights from Biology 

Biological insights for the model behavior can be ob- 
tained by drawing comparisons between the multi-layer 
information system and animal herding behavior [69H7T] . 
The combination of the averaging update rule and the de- 
cision threshold forms the mathematical framework for 
herding or social conformity [72j [73] in the sense that an 
animal can act based on inferences from the information 
of other animals. In humans, models of such information 
diffusion processes are built on a long history of empiri- 
cal work in sociology known as diffusion of innovations 
[351 [7! [IS]. 

For a wide range of animals including humans, group 
decisions to move (for food or travel) often depend on so- 
cial interactions among group members [69], only a few 
of whom have pertinent information (e.g, food location 
or migration route). So-called 'informed' individuals cor- 
respond to our global broadcast, while 'uninformed' in- 
dividuals correspond to our agents on the social network. 
Unlike work investigating individuals who have no deci- 
sion preference [71] , our update rule hard-codes the fact 
that agents in our system have a preference for retaining 
decision states. Such an inertia is consistent with the 
observation that individual beliefs are continually evolv- 
ing variables that depends both on past beliefs [76] and 
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FIG. 7. Facilitation and Hinderance by the Social Network in Case 2. (A) The difference between the average 
information state at d=0.0001 ({Sj (£, 0.0001)}) and that at four other increasing values of the influence parameter showing 
predominantly hinderance (d = 1, top; d = 0.1 second row), a combination of hinderance and facilitation (d = 0.01, third 
row), and predominantly facilitation (d = 0.001, fourth row) by the social network exchange. (B) The larger surface from 
which the four curves in A are drawn. Here we compare (Sj(t, 0.0001)) to d ranging in 0.0001 < d < 1. Color indicates 
(^(£,0.0001)) -(^(M)). 



newly acquired information, and in particular become 
less malleable as time passes [771 [ZE] • 

Insights from Statistical Mechanics 

Physical insights for the model behavior may be ob- 
tained by drawing an analogy with the non-equilibrium 
statistical mechanics of a spin system [79j|80]. Our global 
broadcast plays the role of a uniform external field, while 
the information and adoption state variables of each 
agent can be thought of as continuous or binary spins 
on a directed lattice [45 . The initial state corresponds 
to all spins being initialized at zero, with the external 
field fixed at unity The interaction between individual 
spins and the field is sampled stochastically, leading to 
a noisy dynamical transition from inaction to action (a 
global attractor). 

Inclusion of the social network corresponds to pairwise, 
directed interactions between spins on a random lattice 
[8TH83] that compete with the external field. Because 
the spins are initialized at zero (opposite the field), the 
social network initially tends to hinder action adoption, 
and in some cases prevents action adoption entirely Here 
inclusion of the network homogenizes the collective be- 
havior because the interaction between spins described 
by the update rule is intrinsically stabilizing, damping 
the opinions of outliers back toward the collective. 

A familiar characteristic from the statistical physics 
of spins on a lattice that is not observed in our model is 
the separation of agents/spins into spatially localized do- 



mains characterized by action or inaction [84, 85 . Two 
potential contributing factors are the damping effect of 
our update rule and the mixing effect of the random so- 
cial network lattice, both of which inhibit local propaga- 
tion of injected information. Our preliminary investiga- 
tions suggest that the update rule is the larger predic- 
tor of behavior and therefore we expect that within our 
general framework, clustering would occur for different 
broadcast and opinion update rules. Exploiting parallels 
with well-understood systems in non-equilibrium statis- 
tical physics, operations research, and graph theory are 
likely to provide pathways for systematically unraveling 
the role of underlying network structure, communication, 
and influence of collective behavior of populations. 

Future Directions 

Update Rule Our focus here is on the collective im- 
pact of individual decision-making for when to evacuate, 
rather than the specifics of transportation and routing. 
We implemented an unbiased rule for opinion updates, in 
which the weight assigned to state variables is indepen- 
dent of the state value and the time since the last update. 
This deliberately avoids destabilizing mechanisms pro- 
moting microscopic propagation that arise in other con- 
texts, such as the spread of infectious diseases [86] . In the 
case of evacuations, it is plausible that recent informa- 
tion may be weighted more heavily or travel preferentially 
along specific paths [S7J [88] in opinion updates, or an 
agent with new information may be more likely to share 
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information on the social network. However, inclusion 
of these effects requires more sophisticated assumptions 
about the individual agents that must be justified with 
cognitive and behavioral data. While extracting influ- 
ence and decision rules from network databases remains 
a challenging problem in model identification, studies in 
behavioral psychology, economics, and risk [89j [90] may 
provide useful insights for more realistic representations 
of how individual opinions are updated and decisions are 
made in the context of social networks. 

Network Topology In this initial investigation, our 
model construct is deliberately chosen to be generic, ab- 
stract, and random, setting a baseline for future work. 
The global broadcast is accessed at random by individ- 
uals, and the topology of the social network is random 
as well. In most situations the social network acts over- 
all as a damping force, homogenizing opinion states and 
delaying action adoption in the community, a behavior 
that has also been identified in experimental studies [91] . 
However, even in this scenario, there are situations in 
which the social network accelerates action adoption. 

In each case, design and optimization could naturally 
play a role in policy decisions for specific scenarios. For 
example, the topology of the social network, both in 
terms of connectivity and rate of information flow, could 
be based on realistic measurements of network traffic as 
measured by cellular communication [92 , Twitter [93] . 
or Facebook [94 . Alternatively, the social network for 
transmission of information pertaining to action adop- 
tion during an evacuation might be chosen to correspond 
to the geospatial layout of neighborhoods in a commu- 
nity [95U98] . For a given topology, global broadcast 
rates and/or transmission to individuals could be tailored 
to the connectivity, and optimized for effectiveness as a 
guideline for policy. The clustering inherent in such mod- 
els will likely decrease the cascade effects seen in our data 

EES E2]. 



Concluding Remarks 

In a broader context, the systematic development of 
a framework for understanding the impact of social net- 
works on collective behavior corresponds to the develop- 
ment of a non-equilibrium statistical mechanics for the 
impact of policy on decision-making in populations. His- 
torically, many fields in the social and life sciences have 
taken a phenomenological approach to estimating the im- 
pact of deliberate policy and other externalities on the 
behavior of populations. In this context, the study of 
social networks corresponds to an underlying statistical 
mechanics for the collective behavior, and must be under- 
stood systematically to obtain predictive behavior of the 
population as well as a characterization of the variability 
within the population. 

It is increasingly recognized, across a broad range of 
fields, that understanding network phenomena is essen- 
tial to characterizing behavior of the system as a whole. 
In the specific context of social systems, interactions be- 
tween individuals can for example give rise to financial 
crashes |101] , political revolutions |102] , successful tech- 
nologies |1Q3] , and cultural market sensations |1Q4] . Con- 
structing a statistical mechanics for these problems is 
only a starting point. Feedback, design, optimization, 
and robustness are all critical ingredients, mandating an 
interdisciplinary approach to developing a reliable, pre- 
dictive framework that is useful for policy. Policy issues 
related to individual decisions including the impact of 
training, identification of leaders, and signatures of stress 
and fatigue, will be important topics of future research. 
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Threshold Decisions Our use of a uniform or random 
threshold rule for action adoption at the individual level 
is a traditional starting point used in the decision-making 
literature [38] to study collective population dynamics. 
Other factors that may influence decision-making may 
incorporate time dependent effect — is an individual more 
likely to adopt an action when the rate of action adoption 
is increasing in the population as a whole? Moreover, 
the coupling between multiple networks - e.g., human 
decisions and transportation systems used for evacuation 
- can play a vital role in the collective behavior [100 (e.g., 
when road congestion prevents the group from evacuating 
effectively). A fundamental question is whether these and 
other feedbacks are a cause or an effect of the underlying 
decision process. 
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